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A new strategy of fitting the coupling constants of the nuclear energy density functional is pro- 
posed, which shifts attention from ground-state bulk to single-particle properties. The latter are 
analyzed in terms of the bare single-particle energies and mass, shape, and spin core-polarization 
effects. Fit of the isoscalar spin-orbit and both isoscalar and isovector tensor coupling constants 
directly to the f^/2 — fiji spin-orbit splittings in ''"Ca, ^®Ni, and **Ca is proposed as a practical 
realization of this new programme. It is shown that this fit requires drastic changes in the isoscalar 
spin-orbit strength and the tensor coupling constants as compared to the commonly accepted values 
but it considerably and systematically improves basic single-particle properties including spin-orbit 
splittings and magic-gap energies. Impact of these changes on nuclear binding energies is also 
discussed. 

PACS numbers: 21.10.Hw, 21.10.Pc, 21.60.Cs, 21.60.Jz 
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I. INTRODUCTION 

In this work we propose and explore two new ideas per- 
taining to the energy density functional (EDF) methods. 
First, we suggest a necessity of shifting attention and fo- 
cus of these methods from ground-state bulk properties 
(e.g. total nuclear masses) to single-particle (s.p.) prop- 
erties, and to look for a spectroscopic-quality EDFs that 
would correctly describe nuclear shell structure. Proper 
positions of s.p. levels are instrumental for good descrip- 
tion of deformation, pairing, particle-core coupling, and 
rotational effects, and many other phenomena. 

On the one hand, careful adjustment of these posi- 
tions were at the heart of tremendous success of phe- 
nomenological mean- field (MF) models, like those of Nils- 
son, Woods-Saxon, or folded Yukawa On the other 
hand, similarly successful phenomenological description 
of nuclear masses, within the so-called microscopic- 
macroscopic method Q, relies on the liquid-drop mass 
surface, which is entirely decoupled from the s.p. struc- 
ture. 

Up two now, methods based on using EDFs, in any 
of its variants like local Skyrme, non-local Gogny, or 
relativistic-mean- field (RMF) [s*] approach, were mostly 
using adjustments to bulk nuclear properties. As a result, 
shell properties were described poorly. After so many 
years of investigations, a further increase in precision and 
predictability of all methods based on the EDFs may re- 
quire extensions beyond forms currently in use. Before 
this can be fully achieved, we propose to first take care of 
the s.p. properties, and come back to precise adjustment 
of bulk properties once these extensions are implemented. 

Second, we propose to look at the s.p. properties of 
nuclei through the magnifying glass of odd-even mass 
differences. This idea has already been put forward in 
a seminal paper by Rutz et al. |^], where calculations 
performed within the RMF approach were presented. In 
this paper we perform analogous analysis for the EDFs 



based on the Skyrme interactions. 

On the one hand, it was recognized long time ago, 
cf., e.g., Refs. 0, @, that the theoretical s.p. energies, 
defined as eigenvalues of the MF Hamiltonian, cannot 
be directly compared to experiment, because they are 
strongly renormalized by the particle-core coupling. On 
the other hand, procedures used to deduce the s.p. ener- 
gies from experiment {1^, ^, ^] require various theoretical 
assumptions, by which these quantities cease to result 
from direct experimental observation. In the past, these 
theoretical and experimental caveats hampered the use of 
s.p. energies for proper adjustments of EDFs. However, 
the odd-even mass differences carry very similar physi- 
cal information to that given by s.p. energies, and have 
advantage of being clearly defined, both experimentally 
and theoretically. 

Indeed, experimental difference in mass between an 
odd nucleus and its lighter even-even neighbor, i.e., the 
particle separation energy, is an easily available and un- 
ambiguous piece of data, which reflects the physical role 
of the s.p. energy, with all polarizations and couplings 
taken into account. Similarly, differences of masses be- 
tween the low-lying excited states in an odd system and 
lighter even-even neighbor may illustrate effective posi- 
tions of higher s.p. states. Physical connections between 
these mass differences and s.p. energies are closest in 
semi-magic nuclei, which will be studied in the present 
paper. 

In theory, the primary goal of the EDF methods is to 
describe ground-state energies of fermion systems, i.e., in 
nuclear-physics applications - masses of nuclei. For odd 
systems, the EDF methods should give masses of ground 
states and of several low-lying excited states; the lat- 
ter being obtained by blocking specific s.p. orbitals. We 
stress here that in an odd system, separate self-consistent 
calculations have to be performed for each of the blocked 
states, so as to allow the system for exploring all possible 
polarizations exerted by the odd particle on the even-even 
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core. Again, this procedure is clearly defined and entirely 
within the scope and spirit of the EDF method, which is 
supposed to provide exact energies of correlated states. 
Note, that although one here calculates total masses of 
odd and even systems, the comparison with experiment 
involves the differences of masses, for which many effects 
cancel out. Therefore, one can confidently attribute cal- 
culated differences of masses to properties of effective s.p. 
energies, with all polarization effects taken into account, 
like in the experiment. Again, in semi-magic nuclei such 
a connection is most reliable. 

Following the above line of reasoning, in the present 
paper we study experimental and theoretical aspects of 
the spin-orbit (SO) splitting between the s.p. states. In 
particular, we analyze the role of the SO and tensor 
MFs in providing the correct values of the SO split- 
tings across the nuclear chart. The paper is organized 
as follows: In Sect. [ll] we briefly recall basic theoretical 
building blocks related to the SO and tensor terms of 
the EDFs and interactions. In Sect. IIIII we discuss in 
details three major sources of the core polarization: the 
mass, shape, and spin polarizabilities. In Sect. IIVI we 
present a novel method that allows for a firm adjustment 
of the SO and tensor coupling constants arguing that 
currently used functionals require major revisions con- 
cerning strengths of both these terms. The analysis is 
based on the /7/2 — /5/2 SO splittings in three key nuclei, 
including spin-saturated isoscalar ^°Ca, spin- unsaturated 
isoscalar ^^Ni, and spin- unsaturated isovector ''^Ca sys- 
tems, and is subsequently verified by systematic calcula- 
tions of the SO splittings in magic nuclei. In Sect. Ivl we 
discuss an impact of these changes on binding energies in 
magic nuclei. Conclusions of the paper are presented in 
Sect. IVll 
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For the time-even, pt, r^, and Jt, and time-odd, Sj, T(, 
and jt, local densities we follow the convention introduced 
in Ref. see also Refs. [ll[ and references cited 

therein. 

In particular, the SO density J is the vector part of 
the spin-current tensor density JJ, i.e., 
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In the context of the present study, the time-even ten- 
sor and SO parts of the EDF, 



■Hso = Co'^'VoV- Jo + C7VV- Ji, 



(8) 
(9) 



are of particular interest. 

In the spherical-symmetry limit, the scalar j'-*'-' and 
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symmetric-tensor Jl^J parts of the spin-current tensor 
vanish, and thus 



Ht = \CiJl{r) + \C(jl{r), 



(10) 



Hso = -C^'Ur)'^-^-CY'J.{r)'^-^, (11) 
ar ar 



II. SPIN-ORBIT AND TENSOR ENERGY 
DENSITIES, MEAN FIELDS, AND 
INTERACTIONS 

We begin by recalling the form of the EDF, which will 
be used in the present study. In the notation defined in 
Ref. [l3| (see Ref. [ll| for more details and extensions), 
the EDF reads 



£ = / d'V7^(r), 



(1) 



where the local energy density 7Y(r) is a sum of the kinetic 
energy, and of the potential-energy isoscalar {t — 0) and 
isovector {t — 1) terms. 



where the SO density has only the radial component, 
Jt = ^Jt{r). Variation of the tensor and SO parts of the 
EDF over the radial SO densities J{r) gives the spherical 
isoscalar {t = 0) and isovector {t — 1) SO MFs, 



l-(^Ci'Mr)-C7'f)L.S. (12) 



which can be easily translated into the neutron {q ~ n) 
and proton (g = p) SO MFs, 
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(Co^-C/)Jo(r)+2Cfj,(r) (13) 



dr 



with 



n{r)^—To + Ho{r)+Hiir), 
Zm 



Ht(r) =7^r™(r) +Wr (r) 



(2) 



(3) 



Although below we perform calculations without assum- 
ing spherical and time-reversal symmetries, here we do 
not repeat general expressions for the SO mean-fields, 
which can be found in Refs. [O, [l3- We also note that, 
in principle, in this general case, one could use different 
coupling constant multiplying each of the three terms in 
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Eq. ([7]). In the present exploratory study, we do not 
implement this possible extension of the EDF, and we 
use unique tensor coupling constants C/, as defined in 
Eq. Q. 

Identical potential-energy terms of the EDF, Eqs. (|4]) 
and (III), are obtained by averaging the Skyrme effective 
interaction within the Skyrme-Hartrcc-Fock (SHE) ap- 
proximation Q. By this procedure, the EDF coupling 
constants Ct can be expressed through the Skyrme-force 
parameters, and one can use parameterizations existing 
in the literature. It is clear that one can study tensor 
and SO effects entirely within the EDF formalism, i.e., 
by considering the corresponding tensor and SO parts of 
the EDF, Eqs. ([8]) and and coupling constants C/ 
and C^"^ , respectively. However, in order to link this 
approach to those based on the Skyrme interactions, we 
recall here expressions based on averaging the zero-range 
tensor and SO forces [H, [13, see also Refs. [ill, E m 
for recent analyses. Namely, in the spherical-symmetry 
limit, one has 

= 1 [teJn{r)Jp{r)+t,{J^{r) - J„{r)Jp{r))] , (14) 
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where in Eq. ([TS]) two different coupling constants, Wq 
and Wi , were introduced following Rcf . [17,] . 
The corresponding SO MFs read 

= ^\l(ito + to)Mr)-{t.^QMr)) (16) 



2Wi 



dr 



L-S, 



[note that in Ref. [l^ , the factor of i was missing at the 
Wo term of Eq. (4)]. By comparing Eqs. ^ and (fT6)) . 
one obtains the following relations between the coupling 
constants: 



C( = 

Og — 
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tc), 



(17) 
(18) 
(19) 
(20) 



For further discussion of the Skyrme forces and their rela- 
tion to tensor components we refer the reader to an exten- 
sive and complete recent discussion presented in Ref. 1^ . 

In this exploratory work, we base our considerations 
on the EDF method and deliberately break the connec- 
tion between the functional ([3]), and the Skyrme central, 
tensor, and SO forces. Nevertheless, in the time-even 
sector, our starting point is the conventional Skyrme- 
force-inspired functional with coupling constants fixed at 
the values characteristic for either SkP 18], SLy4 p^ . 
or SkO [20] Skyrme parameterizations. However, poorly 
known coupling constants in the time-odd sector (those 



which are not related to the time-even ones through the 
local-gauge invariance [To| ) are fixed independently of 
their Skyrme-force values. For this purpose, the spin cou- 
pling constants are readjusted to reproduce empirical 
values of the Landau parameters, according to the pre- 
scription given in Refs. [2l|, [S^], and Q'^" are set equal 
to zero. These variants of the standard functionals are 
below denoted by SkPi, SLy4L, and SkOi. 

Strictly pragmatic reasons, like technical complexity 
and lack of firm experimental benchmarks, made the 
majority of older Skyrme parameterizations simply dis- 
regard the tensor terms, by setting C/ = 0. Recent 
experimental discoveries of new magic shell-openin g's in 
neutron-reach light nuclei, e.g., around iV = 32 [23lT2^. 
and their subsequent interpretation in terms of tensor in- 
teraction within the shell-model jl^, [2^ , caused a re- 
vival of interest in tensor terms within the MF approach 
[H, [H, m, H, which is naturaUy tailored to 

study s.p. levels. Indeed, as shown in Ref. _15j, the ten- 
sor terms mark clear and unique fingerprints in isotonic 
and isotopic evolution of s.p. levels and, in particular, in 
the SO splittings. 

In this paper, we perform systematic study of the SO 
splittings. The goal is to resolve contributions to the SO 
MF dni) due to the tensor and SO parts of the EDF, 
and to readjust the corresponding coupling constants C/ 
and C^"^ ■ It is shown that this goal can be essentially 
achieved by studying the /7/2~/5/2 SO splittings in three 
key nuclei: ^''Ca, ''^Ca, and ^®Ni. Before we present in 
Sec. IIVI details of the fitting procedure and values of the 
obtained coupling constants, first we discuss effects of the 
core polarization and its influence on the calculated and 
SO splittings. 



III. CORE-POLARIZATION EFFECTS 

In spite of the fact that the s.p. levels belong to basic 
building blocks of the MF methods, there is still a vivid 
debate concerning their physical reality. The question 
whether they constitute only a set of auxiliary quanti- 
ties, or represent real physical entities that can be in- 
ferred from experimental data, was never of any con- 
cern for methods based on phenomenological one-body 
potentials. Indeed, these potentials were bluntly fitted 
directly to reproduce the s.p. levels deduced, in one way 
or another, from empirical information around doubly- 
magic nuclei, see, e.g., Refs. d, [S^]- In turn, these 
potentials also appear to properly (satisfactorily) repro- 
duce the one-quasiparticle band-heads in open-shell nu- 
clei, see, e.g., Ref. [33 ■ This success seems to legitimate 
the physical reliability of the theoretical s.p. levels within 
the microscopic-macroscopic approaches. 

The debate concerns mostly the self-consistent MF ap- 
proaches based on the EDF methods or two-body ef- 
fective interactions. The arguments typically brought 
forward in this context underline the fact that the self- 
consistent MFs are most often tailored to reproduce bulk 



4 



nuclear properties like masses, densities, radii, and cer- 
tain properties of nuclear matter. Consequently, the un- 
derlying interactions are non-local with effective masses 
m*/m ~ 0.8 [131, what in turn artificially lowers the den- 
sity of s.p. levels around the Fermi energy. It was pointed 
out by several authors 0, H, [sH that restitution of phys- 
ical density of the s.p. levels around the Fermi energy 
can be achieved only after inclusion of particle-vibration 
coupling, i.e., by going beyond MF. 

Within the effective theories, however, these arguments 
do not seem to be fully convincing. First of all, an ef- 
fective theory with properly chosen data set used to fit 
its parameters should automatically select physical value 
of the effective mass and readjust other parameters (cou- 
pling constants) to this value. Examples of such implicit 
scaling are well known for the SHF theory including: (i) 
explicit effective mass scaling of the coupling constants 
C and through the fit to the spectroscopic Landau 
parameters [2l|; (m) direct dependence of the isovector 
coupling constants and Cf on the isoscalar effective 
mass through the fit to the (observable) symmetry en- 
ergy strength [36j : (iii) numerical indications for the m* 
scaling of the SO interaction inferred from the f 7/2-^3/2 
splittings in ^ ~ 44 nuclei [11] . The differences between 
various parameterizations of the Skyrme-force (or func- 
tional) parameters rather clearly suggest that such an 
implicit m* dependence of functional coupling constants 
is a fact, which is, however, neither well recognized nor 
understood so far. 

Secondly, the use of effective interaction with parame- 
ters fitted at the MF level is not well justified in beyond- 
MF approximation. Rather unavoidable double count- 
ing results in such a case, and quantitative estimates of 
level shifts resulting from such calculations need not be 
very reliable. It is, therefore, quite difficult to accept 
the viewpoint that unsatisfactory spectroscopic proper- 
ties of, in particular, modern Skyrme forces can be cured 
solely by going beyond MF. On the contrary, the magni- 
tude of discrepancies between the SHF and experimental 
s.p. levels (see below) rather clearly suggest that: (?) the 
data sets used to fit the force (or functional) parameters 
are incomplete and (ii) the interaction/functional should 
be extended. This is exactly the task undertaken in the 
present exploratory work. We extend the conventional 
EDFs based on Skyrme interactions by including tensor 
terms, and fit the corresponding coupling constants to 
the SO splittings rather than to masses. The prelimi- 
nary goal is to improve spectroscope properties of func- 
tional, even at the expense of the quality in reproducing 
the binding energies. 

Empirical energy of a given s.p. neutron orbital can be 
deduced from the difference between the ground-state en- 
ergy of the doubly- magic core with N neutrons, Eo{N), 
and energies, Ep{N + 1) or Eh{N — 1), of its odd neigh- 
boring isotopes having a single-particle (p) or single-hole 



(h) occupying that orbital, i.e., 

ep(iV) = Ep[N + l)-Eo{N), (21) 
eh{N) = Eo{N)-Eh{N -I), (22) 

see, for example, Refs. [3, H, [1] and referenced cited 
therein. Note that the total energies above are nega- 
tive numbers and decrease with increasing numbers of 
particles, Eh{N - 1) > Efi{N) > Ep{N + 1). Similarly, 
we can define a measure of the neutron shell gap as the 
difference between the lowest single-particle and highest 
single-hole energy, 

Aegap(iV) = ep{N)-eh{N) (23) 
= EpiN + 1) + EhiN - 1) - 2EoiN). 

Single-particle energies of proton orbitals and proton 
shell gaps are defined in an analogous way. 

Consistently with the empirical definitions, in the 
present paper, the same procedure is used on the the- 
oretical level [j]. It means that we determine the total 
energies of doubly-magic cores and their odd neighbors 
by using the EDF method, and then we calculate the s.p. 
energies as the corresponding differences (PT|) or (22). In 
this way, we avoid all the ambiguities related to questions 
of what the s.p. energies really are and how they can be 
extracted from data. Our methodology simply amounts 
to a specific way of comparing measured and calculated 
masses of nuclei. In odd nuclei, one has to ensure that 
particular s.p. orbitals are occupied by odd particles, but 
still the calculated total energies correspond to masses of 
their ground and low-lying excited states. Since here we 
consider only doubly-magic nuclei and their odd neigh- 
bors, the pairing correlations are neglected. 

In our calculations, both time-even and time-odd po- 
larizations exerted by an odd particle/hole on the doubly- 
magic core are evaluated self-consistently. In order to 
discuss these polarizations, let us recall that the extreme 
s.p. model, in which ground-state energies are sums of 
fixed s.p. energies of occupied orbitals, is the model with 
no polarization of any kind included. In this model, the 
differences of ground-state energies, Eqs. (PT|) and ([^ . 
are, of course, exactly equal to the model s.p. energies. 
Having this background model in mind, wc can distin- 
guish three kinds of polarization effects: 

A. Mass polarization effect (time-even) 

This effect can be understood as the self-consistent 
rearrangement of all nucleons, which is induced by an 
added or subtracted odd particle, while the shape is con- 
strained to the spherical point and the time-odd mean 
fields are neglected. One has to keep in mind, that 
none of the s.p. orbitals of a spherical multiplet is spher- 
ically symmetric, and therefore, when a particle in such 
a state is added to a closed core, the resulting state of an 
odd nucleus cannot be spherically symmetric, and thus 
cannot be self-consistent. A self-consistent solution can 
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only be realized within the so-called filling approxima- 
tion, whereupon one assumes the occupation probabili- 
ties of all orbitals belonging to the spherical multiplet of 
angular momentum j and degeneracy 2j + I, to be equal 
tov] = l/{2j + l). 

The mass polarization effect is exactly zero when eval- 
uated within the first-order approximation. Indeed, as 
expressed by the Koopmans theorem [37] , the linear term 
in variation of the total energy with respect to adding or 
subtracting a particle is exactly equal to the bare single- 
particle energy of the core. However, the obtained results 
do not agree with Koopmans theorem. This is illustrated 
in Fig.[Tl where we compare bare and polarized s.p. ener- 
gies of the i^f 7/2 and 1^^3/2 orbitals in ■*°Ca. One can see 
that for both orbitals there is a quite large and positive 
mass polarization effect. 

The reason for this disagreement lies in the fact that, 
because of the center-of-mass correction, the standard 
EDF calculations are not really variational with respect 
to adding or subtracting a particle. Indeed, in these cal- 
culations, the total energy of an A-particle system is cor- 
rected by the center-of-mass correction [38,, ,39,] , 
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(24) 



where i^kin is the average kinetic energy of the system. 
The factor of -j is not varied when constructing the stan- 
dard mean-field Hamiltonian, and therefore, the Koop- 
mans theorem is violated. Had this variation been in- 
cluded, it would have shifted the mean-field potential and 
thus all bare s.p. energies in given nucleus by a constant, 



■^E^in- 



(25) 



In ''°Ca, this shift equals to 0.40 MeV and almost ex- 
actly corresponds to the entire mass polarization effect 
shown in Fig. [T] The remaining mass polarization effect 
can be attributed to higher-order polarization effects (be- 
yond linear) and to the non-selfconsistency of the filling 
approximation. 

A few remarks about the shift in Eq. (|25p are here in 
order. First, the effect is independent on whether the 
two-body or one-body center-of-mass correction [s^ is 
used, and on whether this is done before or after vari- 
ation. Only the detailed value of the shift may depend 
on a particular implementation of the center-of-mass cor- 
rection. Second, the shift induces an awkward result of 
the mean-field potential going asymptotically to a posi- 
tive constant and not to zero. Although this may seem 
to be a trivial artifact, which does not influence the s.p. 
wave functions and observables, it shows that the stan- 
dard center-of-mass correction should be regarded as an 
ill-defined theoretical construct. This fact shows up as an 
acute problem in fission calculations [40| . Third, when- 
ever the bare s.p. energies are compared to empirical 
data, this shift must by taken into account. Alternatively, 
as advocated in the present study, one should compare 
directly the calculated and measured mass differences. 
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FIG. 1: Bare, mass polarized, mass and shape polarized 
(time-even) and mass, shape and spin polarized s.p. ener- 
gies of the uf7/2 particle and the fd^yj i^ole in ''"Ca. The 
polarized s.p. levels were deduced from binding energies in 
one particle/hole nuclei according to the formulas ((21} and 
(|22p . Note that attractive (repulsive) polarization corrections 
to the binding energies shift the particle state down (up) and 
the hole state up (down), respectively. 



Finally, one should note that the shift is irrelevant when 
differences of the s.p. energies are considered, such as the 
SO sphttings discussed in Sec. IIVI 



B. Shape polarization effect (time-even) 

This effect is well known both in the MF and particle- 
vibration-coupling approaches. Within a deformed MF 
theory (with the time-odd mean fields neglected), it cor- 
responds to a simple fact that the s.p. energies (eigen- 
energies in a deformed potential) depend on the defor- 
mation in specific way, which is visualized by the stan- 
dard Nilsson diagram |lj] . Indeed, in an axially deformed 
potential, a spherical multiplet of angular momentum 
j splits into j + 1/2 orbitals according to moduli of 
the angular-momentum projections K — \mj\. Unless 
j = 1/2, for prolate and oblate deformations, orbitals 
with K — 1/2 decrease and increase in energy, respec- 
tively, while those with maximum K — j behave in an 
opposite way. Therefore, both for prolate and oblate de- 
formations, and for j > 1/2, the lowest orbitals have the 
energies that are lower than those at the spherical point. 
Hence, a. j > 1/2 particle added to a doubly-magic core 
always polarizes the core in such a way that the total 
energy decreases. On the other hand, the energy of a 
j = 1/2, K = 1/2 orbital does not depend on deforma- 
tion (in the first order), and thus such an orbital does 
not exert any shape polarization (in this order). 

Exactly the same result is obtained in a particle- 
vibration-coupling model, in which a j > 1/2 particle can 
be coupled with either 0+ or 2+ state of the core, [j®0^]j 
or (g)2"'"]j, and the repulsion of these two configurations 
decreases the energy of the ground state with respect to 
the unperturbed spherical configuration [j (E) As 
before, for j = 1/2, the configuration [j (g) 2+]^ does not 
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exist, and the ground state is not lowered. 

The above reasoning can be repeated for hole states, 
with the result that the j > 1/2 holes added to the 
doubly-magic core always polarize the core in such a way 
that the total energy also decreases. As a consequence, 
the shape polarization effect decreases the s.p. energies 
of particle states (|21l) and increases those of hole states 
(P^ . and thus decreases the shell gap ([M]) . This effect is 
clearly illustrated in Fig. [TJ 

C. Spin polarization effect (time-odd) 

When an odd particle or hole is added to the core, and 
the time-odd fields are taken into account, it exerts po- 
larization effects both in the time-even (mass and shape) 
and time-odd (spin) channels. It means that a non-zero 
average spin value of the odd particle induces a time-odd 
component of the mean field, which influences average 
spin values of all particles, leading to a self-consistent 
amplification of the spin polarization. 

It should be noted at this point that the spin polar- 
ization effect dramatically depends on the assumed sym- 
metries and choices made for the occupied orbital, see 
discussion in the Appendix of Ref. [4l|. Indeed, without 
the time-odd fields, in order to occupy the odd parti- 
cle or hole, one can use any linear combination of states 
forming the Kramers-degenerate pair. The total energy 
is independent of this choice, because the time-even den- 
sity matrix does not depend on it. This allows for mak- 
ing specific additional assumptions about the conserved 
symmetries, e.g., in the standard case, one assumes that 
the odd state is an eigenstate of the signature (or sim- 
plex) symmetry with respect to the axis perpendicular 
to the symmetry axis. However, in order to fully allow 
for the spin polarization effects through time-odd fields, 
one has to release all such restrictive symmetries and al- 
low for alignment of the spin of the particle along the 
symmetry axis. This requires calculations with broken 
signature symmetry and only the parity symmetry be- 
ing conserved. Therefore, calculations of this kind are 
more difficult than those performed within the standard 
cranking model. 



D. Total polarization effect 

In Tables H] and |TT] we list the neutron s.p. energies 
calculated in six doubly-magic nuclei for the SLy4i in- 
teraction. The bare s.p. energies (a) are compared to 
those calculated from total energies, Eqs. (|2T|) and ((22)l . 
with the mass and shape (b) or mass, shape, and spin 
(c) polarizations included. In order to remove ambigu- 
ities associated with occupancy of the valence particle 
(hole) , binding energies of odd- A nuclei were calculated 
by blocking the lowest (highest) K = j orbitals at oblate 
(prolate) shape for particle (hole) orbitals. The blocked 
orbitals were selected by performing cranking calculation 



with angular-frequency vector parallel to the symmetry 
axis. Such a cranking does not affect total energy or wave 
function, but splits spherical multiplets into orbitals hav- 
ing good projections of the angular momentum on the 
symmetry axis. Calculations were performed by using 
the code HFODD (v2.30a) [Ulillii,!!! for the spher- 
ical basis of Nsh — 14 harmonic-oscillator shells. 

As seen by comparing columns (b) and (a) of Tables |T] 
andim the energy shifts caused by the time-even polar- 
ization effects with respect to bare s.p. spectra are almost 
always positive, both for particle and hole states. A few 
exceptions occur only for large-j unfavored (j = £ — 1/2) 
SO partners in heavy nuclei. These shifts clearly de- 
crease in magnitude with increasing mass, from about 
IMeV in ^^ O to b elow 0.25 MeV in ^°^Ph. As dis- 
cussed in Sec. IIII Al they are mainly caused by the mass- 
polarization effects related to the center-of-mass correc- 
tion. Indeed, shifts of s.p. energies calculated for 
the six doubly-magic nuclei of Tables HI and HH are 0.87, 
0.40, 0.36, 0.20, 0.14, and 0.09 MeV, respectively. 

It is also clearly visible that shifts of particle states 
are systematically smaller than those of hole states, i.e., 
the time-even polarizations tend to slightly decrease shell 
gaps. 

The time-odd polarization effects systematically shift 
the hole states down and particle states up in energy, i.e., 
they result in an increase of shell gaps, cf. also Fig. [T] 
This result is at variance with that obtained within 
the RMF approach 4], where the time-odd fields corre- 
sponded to magnetic properties driven by the Lorentz in- 
variance, while here they are determined by experimental 
values of the Landau parameters [111, [l^ . We note here 
that in recent derivations of the time-odd coupling con- 
stants within the relativistic point-coupling model [46| . 
one obtains values of the Landau parameters compat- 
ible with experimental values. Shifts of s.p. energies 
due to the time-odd polarization effects also decrease 
with mass, from about — 0.7(-|-0.5) MeV in -'^^O to below 
-0.1(0.15) McV in ^ospb for hole (particle) states. 

The total effect of combined time-even and time-odd 
polarizations results in adding up the shifts for particle 
states and subtracting those for hole states. In this way, 
the total shifts of particle and hole states become mostly 
positive and (apart from light nuclei) comparable in mag- 
nitude, with quite small net effects on shell gaps. They 
also decrease with increasing mass, from up to 1.5 MeV in 
1*^0 to below 0.25 MeV in ^"Spb. Altogether, polarization 
effects turn out to be significantly smaller than those ob- 
tained in previous estimates. Although in quantitative 
analysis they cannot at all be neglected, discrepancies 
with experimental data (last columns in Tables U and |TT| 
are still markedly larger in magnitude. Therefore, bare 
s.p. energies can be safely used, at least in all studies that 
do not achieve any better overall agreement with data. 
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TABLE L Neutron s.p. energies in ^^O, *'^Ca, and **Ca (in MeV). The columns show: (a) bare unpolarized s.p. energies in 
doubly-magic cores, (b) self-consistent s.p. energies obtained from binding energies in one-particle/hole nuclei, Eqs. (|2ip and 
(|22|) . with time-even mass and shape polarizations included, (c) as in (b), but with time-odd spin polarizations included in 
addition, and (d) experimental data taken from Ref. @]. Time-even, (b) — (a), time-odd, (c) — (b), and total, (c) — (a) polarizations 
are also shown, along with the differences between the self-consistent and experimental spectra, (c) — (d). Positive s.p. energies 
are shown only to indicate that particular orbitals are unbound in calculations; their values are only very approximately related 
to positions of resonances. All results have been calculated using the Sly4i, functional. 







bare 


T-even 


T-even 


T-even 


T-odd 


total 


exp. 


theory 










pol. 


& T-odd 


pol. 


pol. 


[8] 


—exp 






(a) 


(b) 


(b)-(a) 


(c) 


(c)-(b) (c)-(a) 


(d) 


(c)-(d) 


IbQ 




-20.57 


-19.61 


0.96 


-20.29 


-0.68 


0.28 


-21.84 


1.55 






-14.54 


-13.55 


0.99 


-13.86 


-0.31 


0.68 


-15.66 


1.80 




li^rf5/2 


-6.75 


-5.83 


0.92 


-5.43 


0.40 


1.32 


-4.22 


-1.21 






-3.78 


-2.79 


0.99 


-2.30 


0.49 


1.48 


-3.35 


1.05 






0.39 


1.19 


0.80 


1.37 


0.18 


0.98 


1.50 


-0.13 


^"Ca 


lz^rf5/2 


-22.01 


-21.55 


0.46 


-21.87 


-0.32 


0.14 


-22.39 


0.52 




2^3-1/2 


-17.25 


-16.93 


0.32 


-17.51 


-0.58 


-0.26 


-18.19 


0.68 






-15.31 


-14.84 


0.47 


-14.98 


-0.14 


0.33 


-15.64 


0.66 




tJ^j7/2 


— y .oo 




U.oD 


Q nn 
— y .uu 


0.22 


0.58 


— O.Dz 


— U.oO 




2^-^3/2 


-5.24 


-4.85 


0.39 


-4.67 


0.18 


0.57 


-6.76 


2.09 




2l^Pl/2 


-3.06 


-2.66 


0.40 


-2.55 


0.11 


0.51 


-4.76 


2.21 




1J^/5/2 


-1.38 


-1.10 


0.28 


-0.99 


0.11 


0.39 


-3.38 


2.39 


«Ca 


1^^4/2 


-22.60 


-22.02 


0.58 


-22.21 


-0.19 


0.39 


-17.31 


-4.90 




2vSxl2 


-17.60 


-17.26 


0.34 


-17.78 


-0.52 


-0.18 


-13.16 


-4.62 




\vd.ii2 


-16.55 


-15.97 


0.58 


-16.02 


-0.05 


0.53 


-12.01 


-4.01 




li^/7/2 


-9.79 


-9.23 


0.56 


-9.34 


-0.11 


0.45 


-9.68 


0.34 




2vp-ii2 


-5.54 


-5.25 


0.29 


-5.12 


0.13 


0.42 


-5.25 


0.13 




2l^Pl/2 


-3.54 


-3.21 


0.33 


-3.11 


0.10 


0.43 


-3.58 


0.47 




l«^/5/2 


-1.33 


-1.25 


0.08 


-1.21 


0.04 


0.12 


-1.67 


0.46 



IV. SPIN-ORBIT SPLITTINGS 

Before proceeding to readjustments of coupling con- 
stants so as to improve the agreement of the SO spUttings 
with data, we analyze the influence of time-even (mass 
and shape) and time-odd (spin) polarization effects on 
the neutron SO splittings. Based on results presented in 
the preceding Section, we calculate the SO splittings as 

Aego = e„fj< - e„£j> . (26) 

Figure [5] shows the SO splittings calculated u sing 
SLy4i — the functional based on the original SLy4 [1^ 
functional with spin fields readjusted to reproduce em- 
pirical Landau parameters according to the prescription 
given in Refs. [U [l^l • Plotted values correspond to re- 
sults presented in Tables |I] and [TTl The results are la- 
beled according to the following convention: open sym- 
bols mark results computed directly from the s.p. spectra 
in doubly-magic nuclei (bare s.p. energies). These bare 
values contain no polarization effect. Gray symbols label 
the SO splittings involving polarization due to the time- 
even mass- and shape-driving effects, i.e., those obtained 
with all time-odd components in the functional set equal 



to zero. Black symbols illustrate fully self-consistent re- 
sults obtained for the complete SLy4i functional. Gray 
and black symbols are shifted slightly to the left-hand 
(right-hand) side with respect to the doubly-magic core 
in order to indicate the hole (particle) character of the 
SO partners. Mixed cases involving the particle-hole SO 
partners are also shifted to the right. 

The impact of polarization effects on the SO splittings 
is indeed very small, particularly for the cases where both 
SO partners are of particle or hole type. Indeed, for these 
cases, the effect only exceptionally exceeds 200 keV, re- 
flecting a cancellation of polarization effects exerted on 
the 2 = £ zt 1/2 partners. The smallness of polariza- 
tion effects hardly allows for any systematic trends to 
be pinned down. Nevertheless, the self-consistent results 
show a weak but relatively clear tendency to slightly en- 
large or diminish the splitting for hole or particle states, 
respectively. 

The situation is clearer when the SO partners are of 
mixed particle (j = £ — 1/2) and hole {j = £ + 1/2) 
character. In these cases, the shape polarization tends 
to diminish the splitting quite systematically by about 
400-500 keV. This behavior follows from naive deformed 
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TABLE II: Same as in TableUbut for *Zr, "^Sn, and 2°«Pb. 







bare 


T-even 


T-even 


T-even 


T-odd 


total 


exp. 


theory 










pol. 


& T-odd 


pol. 


pol. 


[8] 


—exp 






(a) 


(b) 


(b)-(a) 


(c) 


(c)-(b) ( 


c)-(a) 


(d) 


(c)-(d) 


9°Zr 




-23.16 


-22.83 


0.33 


-22.94 


-0.11 


0.22 


-14.76 


-8.18 




11^/5/2 


-17.07 


-16.72 


0.35 


-16.74 


-0.02 


0.33 


-13.05 


-3.69 




2J^P3/2 


-17.52 


-17.30 


0.22 


-17.44 


-0.14 


0.08 


-12.74 


-4.70 




21^^1/2 


-15.46 


-15.26 


0.20 


-15.35 


-0.09 


0.11 


-12.37 


-2.98 




lt'ff9/2 


-12.08 


-11.75 


0.33 


-11.81 


-0.06 


0.27 


-11.69 


-0.12 




2l^d5/2 


-6.73 


-6.59 


0.14 


-6.52 


0.07 


0.21 


-7.20 


0.68 




3j/Si/2 


-4.93 


-4.70 


0.23 


-4.44 


0.26 


0.49 


-5.78 


1.34 




2^^3/2 


-3.99 


-3.62 


0.37 


-3.59 


0.03 


0.40 


-4.77 


1.18 






-3.75 


-3.75 


0.00 


-3.74 


0.01 


0.01 


-4.62 


0.88 


i=*2Sn 


2l/d5/2 


-11.72 


-11.48 


0.24 


-11.56 


-0.08 


0.16 


-9.10 


-2.46 






-9.46 


-9.28 


0.18 


-9.59 


-0.31 - 


-0.13 


-7.55 


-2.04 






-7.66 


-7.30 


0.36 


-7.33 


-0.03 


0.33 


-7.42 


0.09 






-9.11 


-8.91 


0.20 


-8.95 


-0.04 


0.16 


-7.17 


-1.78 






-2.01 


-2.00 


0.01 


-1.95 


0.05 


0.06 


-2.29 


0.34 




3t'P3/2 


0.17 


0.26 


0.09 


0.31 


0.05 


0.14 


-1.31 


1.62 




11^^9/2 


0.95 


0.79 


-0.16 


0.77 


-0.02 - 


-0.18 


-0.91 


1.68 




3^'Pl/2 


0.97 


1.08 


0.11 


1.12 


0.04 


0.15 


-0.72 


1.84 






0.82 


0.84 


0.02 


0.88 


0.04 


0.06 


-0.35 


1.23 


208pj^ 


2iyfr/2 


-12.02 


-11.85 


0.17 


-11.90 


-0.05 


0.12 


-9.96 


-1.94 




l'^«13/2 


-9.52 


-9.29 


0.23 


-9.30 


-0.01 


0.22 


-8.92 


-0.38 




3t'P3/2 


-9.23 


-9.09 


0.14 


-9.17 


-0.08 


0.06 


-8.12 


-1.05 




2v]^l2 


-9.03 


-8.89 


0.14 


-8.91 


-0.02 


0.12 


-7.78 


-1.13 




3j/pi/2 


-8.11 


-8.01 


0.10 


-8.06 


-0.05 


0.05 


-7.72 


-0.34 




2vg<)/2 


-3.19 


-3.19 


0.00 


-3.16 


0.03 


0.03 


-3.73 


0.57 




iviii/2 


-1.53 


-1.65 


-0.12 


-1.67 


-0.02 - 


-0.14 


-3.11 


1.44 




3ud5/2 


-0.50 


-0.46 


0.04 


-0.43 


0.03 


0.07 


-2.22 


1.79 




4l^Si/2 


0.56 


0.65 


0.09 


0.80 


0.15 


0.24 


-1.81 


2.61 




2ug^/2 


0.08 


0.10 


0.02 


0.11 


0.01 


0.03 


-1.35 


1.46 




3l/d3/2 


0.69 


0.76 


0.07 


0.78 


0.02 


0.09 


-1.33 


2.11 



Nilsson model picture where the highest-if members of 
the j = £—1/2 [j = i+1/2) multiplet slopes down (up) as 
a function of the oblate (prolate) deformation parameter. 
As discussed in the previous Section, the time-odd fields 
act in the opposite way, tending to slightly enlarge the 
gap. The net polarization effect does not seem to exceed 
about 300 keV. In these cases, however, we deal with large 
I orbitals having also quite large SO splittings of the 
order of ~8MeV. Hence, the relative corrections due to 
polarization effects do not exceed about 4%, i.e., they 
are relatively small - much smaller than the effects of 
tensor terms discussed below and the discrepancy with 
data, which in Fig. [2] is indicated for the neutron If SO 
splitting in ""^Ca. These results legitimate the direct use 
of bare s.p. spectra in magic cores for further studies 
of the SO splittings, which considerably facilitates the 
calculations. 



As already mentioned, empirical s.p. energies are es- 
sentially deduced from differences between binding en- 
ergies of doubly-magic core and their odd-^ neighbors. 
Different authors, however, use also one-particle transfer 
data, apply phenomenological particle-vibration correc- 
tions and/or treat slightly differently fragmented levels. 
Hence, published compilations of the s.p. energies, and 
in turn the SO splittings, differ slightly from one another 
depending on the assumed strategy. The typical uncer- 
tainties in the empirical SO splittings can be inferred 
from Table IIIIl which summarizes the available data on 
the SO splittings based on three recent s.p. level compi- 
lations published in Refs. 0,[i,[i|- 

Instead of large-scale fit to the data (see, e.g., Refs. [l^, 
Isijl). we propose a simple three-step method to adjust 
three coupling constants Cq^"^, Cq , and C/. The entire 
idea of this procedure is based on the observation that 
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Nucleus 


orbitals 


Ref. [7] Ref. [8] Ref. f9j 


16q 




- ^^1^1/2 


6.17 


6.18 


— 






- V^d-il2 


5.08 


5.72 


5.08 




Trip- 2 


- ^ip7'2 


6.32 


6.32 


— 






- 7rld3/2 


5.00 


4.97 


5.00 






- J^2pi/2 


2.00 


2.00 


1.54 






- J^l/5/2 


4.88 


5.24 


5.64 






- 


6.00 


6.75 


6.75 






- 7r2pi/2 


2.01 


1.72 


1.69 




7rl/7/2 


- 7rl/5/2 


4.95 


5.41 


6.05 






- TTldg/j 


6.00 


5.94 


6.74 






- 2^2^1/2 


— 


1.67 


1.77 






- J^l/5/2 





8.01 


8.80 






- i/ld-/l. 





5.30 


3.08 




7I"2p3/2 


- 7I-2pi/2 





2.14 


1.77 




7rl/7/2 


- 7^1/5/2 





4.92 


— 






- 7rld-/^2 





5.01 


5.29 


56Ni 


1^2^3/2 


- J^2pi/2 


— 


1.88 


1.12 




^^1/772 


- J^l/5/2 





6.82 


7.16 




1-2^3/2 


- 'r2pi/2 





1.83 


1.11 




'^l/7~/2 


- 7rl/5/2 





7.01 


7.50 






- 2^2^3/2 


— 


2.43 


— 




^^15^/2 


- Z^lff7/2 


— 


7.07 


— 






- ^^2^7/2 





0.37 


— 




^^1/7^2 


- ^1/7/2 





1.71 


— 




7r2d5/2 


— 7r2ci3/2 





2.03 


— 




71-159/2 


- 7rlp7/2 





5.56 


— 




7r2p-/2 


- ^2^7/2 


— 


1.50 


— 




'^l/7'/2 


- ^1/7/2 





4.56 


— 






- J^2d3/2 


1.93 


— 


1.93 






- 2^157/2 


7.00 


— 


7.00 




^15972 


- 7^1^7/2 


6.82 


— 


6.82 




'^2p3/2 


- ^2p3/2 


2.85 


— 


2.85 




v1jTI2 


- 2^2/5/2 


2.00 


1.94 


— 




1^3^3/2 


- J^3pi/2 


0.81 


0.59 


1.15 




^1^7l/2 


- i^l/l9/2 


6.53 


6.51 


6.68 




^2d-/2 


- 1^2^3/2 


1.65 


1.93 


1.66 




7r2d5/2 


- •^2d3/2 


1.48 


1.83 


1.75 




^13972 


- 'rlS'7/2 


6.08 


5.33 


6.08 



^Pb 



1^3^5/2 - 

v1gsji2 - 
'^2/,72 - 

7r2/7/2 - 
7r3p3/2 - 

vrl/i-^j 
^2rf7/2 - 



v'id-j,/2 

■ J^2g7/2 
I^llll/2 

- vZPy2 

■ ^2/57^2 

- 7r2/5/2 

- 'r3pi/2 

- 7rl/l9/2 

- ^2^3/2 



0.97 
2.50 
5.84 
0.90 
2.13 
1.93 
0.85 
5.56 
1.68 



0.89 
2.38 
5.81 
0.90 
2.18 
2.02 
0.45 
5.03 
1.62 



0.97 
2.50 
6.08 
0.89 
1.87 
1.93 
0.52 
5.56 
1.46 



TABLE III: Empirical SO splittings. Compilation is based on 
the empirical s.p. levels taken from Refs. liiSj, i9j. 
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FIG. 2: Neutron SO splittings (|26p calculated using the 
SLy4_L functional. White, gray, and black symbols mark bare, 
mass and shape polarized (time-even), and mass, shape, and 
spin polarized (time-even and time-odd) results, respectively. 
Results for hole (particle and particle-hole) orbitals are shifted 
to the left (right) with respect to the core (open symbols) . A 
typical discrepancy with experiment is shown by the arrow in 
«Ca. 



the empirical I/7/2 — I/5/2 SO splittings in ^°Ca, ^^Ni, 
and ^^Ca form very distinct pattern, which cannot be 
reproduced by using solely the conventional SO interac- 
tion. 

The readjustment is done in the following way. First, 
experimental data in the spin-saturated (SS) nucleus 
^•^Ca are used in order to fit the isoscalar SO coupling 
constant Cj'^ ■ One should note that in this nucleus, the 
SO splitting depends only on C^"^ , and not on Cq (be- 
cause of the spin saturation), nor on C^"' (because of 
the isospin invariance at N = Z), nor on (because of 
both reasons above). Therefore, here one experimental 
number determines one particular coupling constants. 

Second, once Cg'"' is fixed, the spin-unsaturated (SUS) 
N = Z nucleus ^^Ni is used to establish the isoscalar 
tensor coupling constant Cq . Again here, because of 
the isospin invariance, the SO splitting is independent 
of either of the two isovector coupling constants, C^"' 
or Ci . Finally, in the third step, ^^Ca is used to adjust 
the isovector tensor coupling constant Cj^. Such a pro- 
cedure exemplifies the focus of fit on the s.p. properties, 
as discussed in the Introduction. 

It turns out that current experimental data, and in 
particular lack of information in ^*Ni, do not allow for 
adjusting the fourth coupling constant, C^"'- For this 
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FIG. 3: Figure illustrates the three-step procedure used to 
fit the isoscalar SO coupling constant C^'' in ''"Ca (up- 
per panel), the isoscalar tensor strength Ci in ^^Ni (middle 
panel), and the isovector tensor strength C( in **Ca (lowest 
panel). These particular calculations have been done for the 
SkP functional, but the pattern is common for all the an- 
alyzed parameterizations including SLy4 and SkO. See text 
for further details. 



reason, in the present study we fix it by keeping the ratio 
of Cj^ /CY'' equal to that of the given standard Skyrme 
force. In the process of fitting, all the remaining time- 
even coupling constants Ct are kept unchanged. Variants 
of the standard functionals obtained in this way are be- 
low denoted by SkPy, SLy4T, and SkOr- When the 
time-odd channels, modified so as to reproduce the Lan- 
dau parameters, are active, we also use notation SkP^T, 
SLy4iT7 and SkO^T- 

For the SkP functional, the procedure is illustrated in 
Fig. [21 We start with the isoscalar N = Z nucleus '*°Ca. 
The evolution of the SO splittings in function of f^'\ 
which is the factor scaling the original SkP coupling con- 
stant Cg^"^, is shown in the upper panel of Fig. [3l As it 
is clearly seen from the Figure, fair agreement with data 
requires about 20% reduction in the conventional SO in- 
teraction strength Cq"^. It should be noted also that the 
reduction in the SO interaction considerably improves 
the 1/7/2 - 1(^3/2 and I/7/2 - 2^3/2 splittings but slightly 
spoils the 2p3/2 — 2pi/2 SO splitting. Qualitatively, sim- 
ilar results were obtained for the SLy4 and SkO interac- 
tions. Reasonable agreement to the data requires ^20% 
reduction of the original C^"^ in case of the SkO inter- 
action and quite drastic ^^35% reduction of the original 
Cj'' in case of the SLy4 force. 

Having fixed C^"^ in ^"Ca we move to the isoscalar nu- 



Skyrme 










force 


[MeVfm^] 




[MeV fm^] 


[MeVfm^] 


SkPr 


-60.0 


3 


-38.6 


-61.7 


SLy4T 


-60.0 


3 


-45.0 


-60.0 


SkOT 


-61.8 


-0.78 


-33.1 


-91.6 



TABLE IV: Spin-orbit C^'' and tensor isoscalar Cq and 
isovector Ci functional coupling constants adopted in this 
work and subsequently used in Figs.|4l[5l and[6l where global 
calculations of the SO splittings are presented. 



cleus ^^Ni. This nucleus is spin- unsaturated and there- 
fore is very sensitive to the isoscalar Cfj tensor coupling 
constant. The evolution of theoretical s.p. levels versus 
Cq is illustrated in the middle panel of Fig. [3l As shown 
in the Figure, reasonable agreement between the empiri- 
cal and theoretical I/7/2 — I/5/2 SO splitting is achieved 
for Cq ^ — 40MeVfm^, which by a factor of about five 
exceeds the original SkP value for this coupling constant. 
It is striking that a similar value of Cq is obtained in the 
analogical analysis performed for the SLy4 interaction. 

Finally, the isovector tensor coupling constant CY'^ is 
established in Nj^Z nucleus ^®Ca. The evolution of the- 
oretical neutron s.p. levels versus is illustrated in the 
lowest panel of Fig. [31 As shown in the Figure, the value 
of CY'^ ~ — 70MeVfm^ is needed to reach reasonable 
agreement for the 11^/7/2 ~ SO splitting in this 

case. For this value of the C/ strength one obtains also 
good agreement for the proton — SO split- 

ting (see Figs. [Hand [51 below), without any further read- 
justment of the CY'' strength. Again, very similar value 
for the C( strength is deduced for the SLy4 force. Note 
also the improvement in the I/7 /2 ~ splitting caused 
by the isovector tensor interaction. Dotted lines show re- 
sults obtained from the mass differences, i.e., with all the 
polarization effects included. 

During the fitting procedure all the remaining func- 
tional coupling constants were kept fixed at their Skyrme 
values. The ratio of the isoscalar to the isovector cou- 
pling constant in the SO interaction channel was locked 
to its standard Skyrme value of Cj'^ /CY"^ = 3. Since 
no clear indication for relaxing this condition is seen (see 
also Figs. [H and \5\ below), we have decided to investi- 
gate the isovector degree of freedom in the SO interac- 
tion (see [l3|) by performing our three-step fitting process 
also for the generalized Skyrme interaction SkO [20|, for 
which CY''/CY'^ ~ -0.78. 

All the adopted functional coupling constants result- 
ing from our calculations are collected in Table HVl Note, 
that the procedure leads to essentially identical SO inter- 
action strengths CY'' for all three forces irrespective of 
their intrinsic differences, for example in effective masses. 
The tensor coupling constants in both the SkP and the 
SLy4 functionals are also very similar. In the SkO case, 
one observes rather clear enhancement in the isovector 
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tensor coupling constant which becomes more negative 
to, most likely, counterbalance the non-standard positive 
strength in the isovector SO channel. 

The functionals were modified using only three specific 
pieces of data on the neutron I/7/2 — I/5/2 SO splittings. 
In order to verify the reliability of the modifications, we 
have performed systematic calculations of the experimen- 
tally accessible SO spHttings. The results are depicted in 
Figs. H El and El for the SkP, SLy4, and SkO function- 
als, respectively. Additionally, Fig. [7] shows neutron and 
proton magic gaps ([24| calculated using the SkP func- 
tional. In all these Figures, estimates taken from Ref. Q 
are used as reference empirical data. 

This global set of the results can be summarized as 
follows: 



Nucleus SLy4 SLy4T SLy4Tmin 

*°Ca -2.197 -1.830 -5.775 

**Ca -1.912 5.039 -0.279 

0.625 15.138 9.127 

^°Zr -1.845 7.492 -3.032 

^^^Sn -0.660 19.898 2.222 

^°^Pb 0.822 24.910 -3.048 

TABLE V: Differences between calculated and experimental 
ground-state energies, -Ecaic — i?cxp, (in MeV) for a set of 
spherical nuclei. Column denoted as SLy4Tmin shows results 
obtained after (local) minimization with respect to parame- 
ters ti,Xi, i = 0,1, 2, 3. See text for details. 



• The £ = 1, P3/2 — Pi/2, SO splittings are slightly 
better reproduced with original rather than modi- 
fied functionals. 

• The Id SO sphttings {^^O, "0.48^^) are rather 
poorly reproduced by both the original and modi- 
fied functionals. These splitting are also subject to 
relatively big empirical uncertainties as shown in 
Table IIIII and, therefore, need not be very conclu- 
sive. In particular, the Id SO splittings in N=Z 
^^O and ^"Ca nuclei deduced from Ref. Q and de- 
picted in the figures show surprisingly large isospin 
dependence. 

• The 2d and 3d splittings are quite well reproduced 
by both the original and modified functionals with 
slight preference for the modified functional, in par- 
ticular for the SLy4 interaction. 

• All the £ > 3 SO splittings are reproduced consid- 
erably better by the modified functionals. 

• Magic gaps are also better (although not fully sat- 
isfactorily) reproduced by the modified functionals. 

Without any doubt the SO splittings are better de- 
scribed by the modified functionals. It should be stressed 
that the improvements were reached using only three ad- 
ditional data points without any further optimization. 
The tensor coupling constants deduced in this work and 
collected in table IIVI should be therefore considered as 
reference values. Indeed, direct calculations show that 
variations in C/ within ±10% affect the calculated SO 
splittings only very weakly. The price paid for the im- 
provements concerns mostly the binding energies, which 
for the nuclei ^^Ni, ^^^Sn, and ^"'^Pb become worse as 
compared to the original values. This issue is addressed 
in the next section. 



V. BINDING ENERGIES 

Parameters of the Skyrme functional have been fitted 
to reproduce several physical quantities, with emphasis 



on the masses of magic nuclei. Therefore it is not sur- 
prising that dramatic modifications of the SO and tensor 
terms of the functional, described in Sec. IIIII while im- 
proving the agreement between the calculated and mea- 
sured single-particle properties, can destroy the quality of 
the mass fit. Hence, it is interesting to know whether this 
disagreement is significant and whether it can be healed 
by refitting the remaining parameters of the functional. 

Table |V] shows differences between calculated and ex- 
perimental (Ref. [131) ground-state energies, i?caic — £^exp, 
(in MeV) for a set of spherical nuclei. Negative val- 
ues mean that nuclei are overbound. Results given in 
the second column, denoted as SLy4, correspond to the 
standard SLy4 [l^ parametrization. The third column, 
denoted as SLy4T, illustrates significant deterioration of 
the quality of fit when parameters Cj^'^, Cq and C/ are 
modified (see Sec. IIIII) . Values presented in the last col- 
umn, SLy4Tmin, wcrc obtained by minimizing the rms of 
relative discrepancies between the calculated and mea- 
sured masses Values of Cj'^ , Cq and C/ were kept fixed 
at their SLy4T values, while minimization was performed 
by varying the remaining parameters of the functional ti 
and Xi. Note that for the standard SLy4 functional, the 
tensor coupling constants are set equal to zero indepen- 
dently of the values of ti and Xi. For the minimization, 
we have used the same methodology, namely, the influ- 
ence of parameters ti and Xi on tensor coupling constants 
was disregarded. 

It should be emphasized that no attempt has been 
made to find the global minimum — the minimization 
was purely local, in the vicinity of the standard SLy4 
values of the parameters ti and Xi. One can see that 
even this very limited procedure can lead to significant 
reduction of discrepancies, down to quite reasonable val- 
ues (with an exception of ^^Ni nucleus). 

It is worth noting that the resulting modifications of 
the ti and Xi parameters turned out to be very small. 
Table IVII shows the values of the ti and Xi parameters 
in the standard SLy4 parametrization (second column) 
and those obtained as the result of the minimization pro- 
cedure (third column). The last column shows relative 
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FIG. 4: Experimental @] (black symbols) and theoretical SO splittings calculated using the original SkP functional (gray 
symbols) and our modified SkPLT functional (white symbols) with the SO and tensor coupling constants given in Table ITVl 
Upper left and right panels show neutron SO splittings for \ow-£ = 1, 2 (p and d) and high-^ > 3 orbitals, respectively. Analogical 
information but for proton SO splittings is depicted in the lower panels. 



changes of parameters (in percent). As one can see, they 
are at most of the order of one percent. Nevertheless, 
even such small changes were sufficient to improve sig- 
nificantly the agreement between calculated and experi- 
mental masses. 

We stress again that the refitting procedure is used 
here only for illustration purposes and the global fit to 
masses must probably include extended functionals and 
improved methodology. For example, the Wigner energy 
correction was not included in the fit, as it was nei- 
ther included in the fit of the SLy4 parametrization. This 
correction alone may change the balance of discrepancies 
obtained for the N = Z and N Z nuclei, and strongly 
impact the results. Systematic studies of these effects 
will be performed in the near future. 



par am. 


SLy4 


SLy4Tmin chauj 


;e (%) 


to 


-2488.913 


-2490.00300 


0.04 


ti 


486.818 


486.78460 


-0.01 


t2 


-546.395 


-545.35849 


-0.19 


t-i 


13777.000 


13767.77776 


-0.07 


Xo 


0.834 


0.83257 


-0.17 


Xl 


-0.344 


-0.34227 


-0.50 


X2 


-1.000 


-0.99798 


-0.20 


X3 


1.354 


1.36128 


0.54 



TABLE VI: Skyrme force parameters ti,Xi, i — 0,1,2,3 of 
the standard SLy4 parametrization (second column) com- 
pared with those obtained from the minimization procedure 
described in the text (third column). The last column shows 
relative change of parameters (in percent). 



VI. CONCLUSIONS 



Applicability of functional-based self-consistent mean- 
field or energy-density-functional methods to nuclear 
structure is hampered by their unsatisfactory s.p. proper- 
ties. This fact seems to be a mere consequence of strate- 
gies used to select datasets that were applied in the pro- 
cess of adjusting free parameters of these effective theo- 
ries. In spite of the fact that the s.p. energies are at the 
heart of these methods, the datasets are heavily oriented 
towards reproducing bulk nuclear properties in large- 
limit, with only a marginal influence of the s.p. levels or 



level splittings in finite nuclei. 

In this work we suggest a necessity of shifting attention 
from bulk to s.p. properties and to look for spectroscopic- 
quality EDF, even at the expense of deteriorating its 
quality in reproducing binding energies. Such a strategy 
requires well-defined empirical input related to the s.p. 
energies, to be used directly in the fitting process. We ar- 
gue that odd-even mass differences around magic nuclei 
not only provide unambiguous direct information about 
nuclear s.p. energies but are also well anchored within 
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FIG. 5: Same as in Fig. [4] but for the SLy4 functional. 
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FIG. 6; Same as in Fig. |4] but for the SkO functional. 



the spirit of the EDF formahsm. Indeed, the theorems 
due to Hochcnberg and Kohn [4§| and Levy [13] , see also 
Refs. 51], imply existence of universal EDF capable, at 
least in principle, treating ground-states energies of nu- 
clei exactly. One can, therefore, argue that this implies 
essentially exact treatment of at least the lowest s.p. lev- 
els forming ground states in one-particle (one-hole) odd- 



A nuclei with respect to even-even cores or, alternatively, 
almost exact description of core-polarization phenomena 
caused by odd single-particle (single-hole). 

An attempt to refit the EDF is preceded by a system- 
atic analysis of the s.p. energies and self-consistent core- 
polarization effects within the state-of-the-art Skyrme- 
force-inspired EDF. Three major sources of core- 
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8 20 28 50 82 126 8 20 28 40 50 82 

FIG. 7; Experimental Q (dots) and theoretical values of magic gaps (|24p . calculated using the original SkP functional (open 
triangles) and the SkPr functional (full triangles) with the SO and the tensor coupling constants from Table. IIVI The gaps 
were computed using the bare unpolarized s.p. spectra. 



polarization, including mass, shape and spin (time-odd) 
effects, are identified and discussed. The analysis is per- 
formed for even-even doubly-magic cores and the lowest 
s.p. states in odd-^ one-particle (hole) nuclei. The discus- 
sion is supplemented by analysis of the s.p. SO splittings. 

New strategy in fitting the EDF is applied to the SO 
and tensor parts of the nuclear EDF. Instead of large- 
scale fit to binding energies we propose simple and in- 
tuitive three-step procedure that can be used to fit the 
isoscalar strength of the SO interaction as well as the 
isoscalar and isovector strengths of the tensor interac- 
tion. The entire idea is based on the observation that the 
/7/2 ~ /5/2 SO splittings in spin-saturated isoscalar '*°Ca, 
spin-unsaturated isoscalar ^^Ni, and spin-unsaturated 
isovector ^^Ca form distinct pattern that can neither be 



understood nor reproduced based solely on the conven- 
tional SO interaction. The procedure indicates a clear 
need for major reduction (from ~20% till ~35% depend- 
ing on the parameterization) of the SO strength and for 
strong tensor fields. It is verified that the suggested 
changes lead to systematic improvements of the func- 
tional performance concerning such s.p. properties like 
the SO splittings or magic gaps. It is also shown that de- 
structive impact of these changes on the binding energies 
can be improved, to a large extent, by relatively small 
refinements of the remaining coupling constants. 

This work was supported in part by the Polish Ministry 
of Science and by the Academy of Finland and University 
of Jyvaskyla within the FIDIPRO programme. 
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